ELSEVIER  Journal  of  Power  Sources  109  (2002)  203-213 


www.elsevier.com /locate /jpowsour 


Modeling  viscosity  and  conductivity  of  lithium  salts  in  y-butyrolactone 

A.  Chagnesa’*,  B.  Carrea,  P.  Willmannb,  D.  Lemordanta 

a Laboratoire  de  Physicochimie  des  Interfaces  et  Milieux  Reactionnels  (EA2048),  Universite  de  Tours , 

Faculte  des  Sciences  et  Techniques,  Parc  de  Grandmont,  F  37200  Tours,  France 
bCNES,  18  Avenue  E.  Belin,  F  31055  Toulouse  Cedex,  France 

Received  20  October  2001;  received  in  revised  form  12  January  2002;  accepted  22  January  2002 


Abstract 

Viscosity  and  conductivity  properties  of  Li-salts  (lithium  tetrafluoroborate  (LiBF4),  lithium  hexafluorophosphate  (LiPF6),  lithium 
hexafluoroarsenate  (LiAsF6),  lithium  bis-(trifluoromethylsulfone)-imide  (LiTFSI))  dissolved  in  y-butyrolactone  (BL)  have  been  investigated. 
The  B-  and  ^-coefficients  of  the  Jones-Dole  (JD)  equation  for  the  relative  viscosity  of  concentrated  electrolyte  solutions  (concentration: 
C  =  0.1-1. 5  M):  ;/r  =  1  +  AC1/2  +  BC  +  DC2,  have  been  determined  as  a  function  of  the  temperature.  The  ^-coefficient  is  linked  to  the 
hydrodynamic  volume  of  the  solute  and  remains  constant  within  the  temperature  range  investigated  (25-55  °C).  The  ZTcoefficient,  which 
originates  mainly  from  long-range  coulombic  ion-ion  interactions,  is  a  reciprocal  function  of  the  temperature.  The  variations  of  the  molar 
conductivity  (A)  with  C  follow  the  cube  root  law  A  =  A0  —  S'C 1//3  issued  from  quasi-lattice  theory  of  electrolyte  solutions.  From  the  Walden 
product  W  =  At]  which  does  not  vary  with  C  and  the  JD  equation,  the  bell  shape  of  the  conductivity-concentration  relationship  is  explained 
and  it  is  shown  that  the  concentration  in  salt  at  the  maximum  of  conductivity  is  linked  to  the  D-coefficient.  Raman  spectroscopy  has  been  used 
as  an  additional  tool  to  investigate  ion  pairing  in  BL.  Ions  pairs  have  been  evidenced  for  LiC104  solutions  in  BL  but  not  for  LiPF6.  As  little 
variations  occur  for  the  ions  pairs  dissociation  coefficient  when  the  salt  concentration  is  increased,  the  cube  root  law  remains  valid,  at  least  in 
the  concentration  range  investigated.  ©  2002  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Viscosity  and  conductivity  of  electrolyte  solutions  are 
related  to  transport  processes  and  play  an  important  role  in 
the  field  of  lithium  batteries  as  the  knowledge  of  the  con¬ 
ductivity  is  necessary  for  the  design  of  Li  cells.  Viscosity  and 
conductivity  measurements,  provide  also  useful  insights  on 
ion  solvation  and  association.  At  the  present  time,  the 
mechanisms  of  ion  transport  in  concentrated  organic  elec¬ 
trolytes  is  not  clearly  understood.  The  lack  of  convenient 
theories  is  evident  and  only  empirical  relations  have  been 
proposed.  As  an  example,  the  extended  Jones-Dole  (JD) 
equation  [1]  for  the  relative  viscosity  of  solutions  (t]v)  has 
been  extensively  used  to  describe  the  viscosity  variations  of 
electrolyte  solutions  when  the  concentration  of  the  salt  (C) 
is  varied: 

>1r  =  >l/%  =  l+ACl/2  +BC  +  DC2  (1) 
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In  this  equation,  t]  and  t]0,  are  the  viscosities  of  the  solution 
and  the  pure  solvent  respectively,  and  A,  B  and  D  are  the 
coefficients.  This  first  term,  in  C  ,  on  the  right  hand  side  of 
Eq.  (1),  is  linked  to  the  interaction  of  a  reference  ion  with  its 
ionic  atmosphere  and  may  be  calculated  by  the  Falkenhagen 
theory  [2]  but,  usually,  this  term  vanished  in  organic  solvents 
when  the  concentration  in  salt  is  raised  above  C  ~  0.05  M. 
The  BC  term  is  predominant  at  C  >  0.05  M  and  has  been 
attributed  to  ion-solvent  interactions  as  well  as  to  volume 
effects  [3].  The  third  term  in  C2,  appears  at  the  highest 
concentrations  in  salt  (0.5-2  M  or  more),  i.e.  when  the  mean 
interionic  distance  decreases  and  becomes  of  the  order  of 
magnitude  of  a  few  solvent  molecules  diameters.  It  is  mainly 
related  to  ion-ion  and/or  ion-solvent  interactions  as  it  does 
not  appear  in  the  case  of  weak  electrolytes  [4]. 

Variations  of  the  molar  conductivity  (A)  of  diluted 
electrolytes  solutions  with  salt  concentration  are  usually 
described  by  the  classical  Debye-Hiickel-On sager  (DHO) 
equation:  [5,6] 

A  =  A°  -  SCl/ 2  (2) 

where  A°  is  the  molar  conductivity  at  infinite  dilution  and  S, 
a  calculable  parameter  which  depends  on  both  the  physical 
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properties  of  the  solvent  and  the  nature  of  the  electrolyte. 
Nevertheless,  experimental  results  indicate  that  the  Eq.  (2)  is 
not  followed  [7,8],  excepted  at  very  low  concentration 
(C  <  0.001  M)  and,  even  in  these  conditions,  the  experi¬ 
mental  slope  (.S’exp)  is  often  different  (higher  or  lower)  than 
the  calculated  one  (Scalc)  [8].  Attempts  to  correct  experi¬ 
mental  result  for  solution  viscosity  and  ion  association  are 
not  successful  and  the  discrepancies  between  theory  and 
experience  remain  [9].  Clearly,  in  concentrated  organic 
electrolytes,  the  conductivity  process  is  different  from  the 
model  embodied  in  the  DHO  equation.  Another  approach 
has  been  provided  by  the  use  of  the  quasi-lattice  theory 
[10,11].  This  theory,  adapted  for  conductivity  [12],  has  been 
applied  successfully  to  the  system  y-butyrolactone  (BL)  + 
LiC104  [7].  According  to  this  theory,  the  variations  of  A  with 
C  are  described  by  the  following  equations: 

A  =  A0/  -  S'C1/3  (3) 

where 

S' =  Si  +  S2A0'  (4) 

_  0.293e2M(2Aa)1/3 
1  47tE0£rkbT 

and 

Aw.W'3  (6) 

7 irj 

In  these  equations,  e0  is  the  permittivity  of  free  space,  gr  the 
relative  permittivity  of  the  solvent,  e  the  electronic  charge,  T 
the  temperature,  kh  the  Boltzmann’s  constant,  Na  the  Avo- 
gadro’s  number,  and  M  the  lattice  parameter  (equivalent  to  a 
Madelung  constant).  A0'  (S  cm-1)  is  the  molar  conductivity 
of  the  salt  at  infinite  dilution  given  by  the  lattice  model.  As 
this  model  is  from  different  ground  from  the  Debye  and 
Hiickel  ones,  the  values  Al"  and  A0,  obtained  by  extrapola¬ 
tion,  respectively,  from  the  C1/3  and  C1/2  law,  may  differ.  As 
the  relaxation  term  Si  is  usually  greater  than  the  electro¬ 
phoretic  term  S2,  the  slope  of  the  A  =  /(C1/3)  relation  is 
mainly  dependent  on  the  lattice  parameter  value  M  and  the 
temperature. 

Propylene  carbonate  (PC)  and  BL  are  popular  dipolar 
aprotic  solvents  widely  used  for  storage  energy  or  electro- 
chromic  devices.  For  practical  applications,  the  lower  visc¬ 
osity  of  BL  over  PC  overrides  its  lower  permittivity.  The 
two  oxygen  atoms  of  BL  act  as  bidentate  ligand  toward 
small  cations  like  Li+  [13]  and  prevent  the  formation  of  non¬ 
conducting  contact  ion-pairs.  As  a  result,  the  moderate 
viscosity  of  BL  and  the  limited  ion  aggregation  are  respon¬ 
sible  for  the  high  conductivity  in  liquid  electrolytes  as  well 
as  in  PMMA-based  gels  [14].  Li-salts  having  large  anions 
lead  to  a  reduction  of  ionic  association  and  thus,  to  an 
increase  in  conductivity  and  solubility:  lithium  tetrafluor- 
oborate  (LiBF4),  lithium  hexafluorophosphate  (LiPF6), 
lithium  hexafluoroarsenate  (LiAsF6),  lithium  bis-(trifluoro- 
methylsulfone)-imide  (LiTFSI)  are  mainly  used  in  practice. 


The  purpose  of  the  present  paper  is  to  examine  whether 
the  preceeding  Eqs.  (1),  (3)-(6)  may  be  applied  to  other  salts 
than  LiC104.  For  this  purpose,  viscosity  and  conductivity 
measurements  of  LiBF4  LiAsF6,  LiPF6  and  LiTFSI  solu¬ 
tions  in  BL  have  been  performed.  All  these  lithium  salts  are 
widely  used  in  the  field  of  lithium  batteries  as  organic 
electrolytes.  Viscosity  measurements  have  been  realised 
from  25  to  55  °C  and  in  a  sufficiently  wide  range  of 
concentration  to  determine  accurately  the  B-  and  D-coeffi- 
cients  of  the  JD  equation  as  a  function  of  the  temperature. 
The  validity  of  the  cube  root  law  for  the  variation  of  the 
conductivity  with  concentration  (Eq.  (3))  has  been  tested  for 
these  electrolyte  solutions.  From  the  temperature  depen¬ 
dence  of  the  viscosity  and  the  conductivity,  the  correspond¬ 
ing  activation  parameters  Earj  and  EaA  have  been 
determined  as  a  function  of  the  salt  concentration.  Raman 
spectroscopy  is  a  useful  tool  to  study  short  range  ion-ion 
interactions,  leading  to  the  possible  formation  of  ion-pairs 
and  higher  ionic  aggregates.  Raman  spectroscopy  is  pre¬ 
sently  used  to  investigate  ion  association  in  BL. 

2.  Experimental 

BL  (purity  >  99%)  was  obtained  from  Aldrich  and  dried 
over  molecular  sieves  prior  to  use.  LiC104  (99%),  LiPF6 
(98%),  LiAsF6  (98%),  LiBF4  (99%),  LiTFSI  (99%)  were 
obtained  from  Fluka  and  dried  in  a  vaccum  oven  before 
dissolution  in  the  solvents.  The  water  content  of  the  solu¬ 
tions  was  <50  ppm  as  indicated  by  Karl  Fisher  titration. 

AC  conductivity  measurements  were  carried  out  using  a 
conductivity  cell  with  platinum  electrodes  and  a  Philips 


Table  1 

Variations  of  the  density  (g  cm~3)  of  LiBF4,  LiPF6,  LiAsF6,  and  LiTFSI 
solutions  in  BL  with  temperature  and  salt  concentration 


LiBF4  (°C) 

0.1  M 

0.2  M 

0.5  M 

0.75  M 

1  M 

1.5  M 

25 

1.130 

1.134 

1.149 

1.161 

1.173 

1.196 

35 

1.120 

1.125 

1.140 

1.152 

1.164 

1.188 

45 

1.112 

1.118 

1.133 

1.144 

1.157 

1.181 

55 

1.105 

1.109 

1.125 

1.136 

1.150 

1.174 

LiPF6  (°C) 

0.2  M 

0.5  M 

0.75  M 

1M 

1.5  M 

25 

1.135 

1.163 

1.186 

1.207 

1.261 

35 

1.128 

1.154 

1.178 

1.198 

1.252 

45 

1.121 

1.146 

1.170 

1.191 

1.244 

55 

1.113 

1.139 

1.162 

1.185 

1.237 

LiAsF6  (°C) 

0.125  M  0.25  M 

0.5  M 

1.5  M 

25 

1.169 

1.157 

1.185 

1.316 

35 

1.131 

1.148 

1.177 

1.308 

45 

1.122 

1.140 

1.169 

1.299 

55 

1.114 

1.133 

1.162 

1.289 

LiTFSI  (°C) 

0.1  M 

0.2  M 

0.5  M 

0.75  M 

1  M 

1.5  M 

25 

1.139 

1.147 

1.183 

1.220 

1.253 

1.318 

35 

1.128 

1.137 

1.173 

1.203 

1.243 

1.303 

45 

1.120 

1.129 

1.157 

1.201 

1.235 

1.297 

55 

1.113 

1.122 

1.147 

1.194 

1.228 

1.285 
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Fig.  1.  Raman  spectrum  of  pure  BL. 


PM6303  impedance  analyser  operating  at  1000  Hz.  Viscos¬ 
ity  measurements  were  obtained  using  an  Ubbelhode  capil¬ 
lary  tube  immersed  in  a  thermostated  bath  (±002  °C)  and  a 
Schott  viscometer(AUS  310).  The  densities  of  solutions, 
required  to  calculate  the  dynamic  viscosities,  were  deter¬ 
mined  using  a  vibrating  tube  densimeter  (Picker)  and  results 
are  reported  in  Table  1. 

The  Raman  spectra  were  recorded  with  a  Bruker  RFS100 
spectrometer.  The  near-infrared  laser  was  operating  at  a 
wavelength  of  1064  nm  with  a  power  on  the  sample  of 
100  mW  and  a  light  spot  of  0.1mm  in  diameter.  The 
scattered  light  was  collected  at  180°  from  the  direction  of 
the  incident  light.  Each  Raman  spectrum  was  the  result  of 
100  scans.  The  germanium  detector  was  cooled  in  liquid 
nitrogen  and  an  OS/2  operating  system  was  used  to  collect 
and  process  data.  Optical  glass  test  tubes  were  used  for 
Raman  measurements.  In  Fig.  1,  is  reported  the  Raman 
spectra  of  pure  BL.  All  the  spectra  have  been  normalised 
using  the  C-H  elongation  of  BL  at  3000  cm  1  because  this 
band  present  a  high  intensity  and  does  not  interfere  with 
others  studied  bands.  In  the  following,  the  contribution  of 
BL  to  Raman  spectra  has  been  removed  by  subtraction. 

3.  Results  and  discussion 

3.1.  Viscosity 

The  variations  of  the  viscosity  at  298  K  of  LiBF4,  LiPF6, 
LiAsF6,  and  LiTFSI  solutions  in  BL  are  plotted  against  the 
salt  concentration  in  Fig.  2.  The  variations  in  viscosity  of 
the  solutions  are  very  similar,  excepted  for  LiBF4  which  is 


the  salt  with  the  smallest  anion  and  the  most  susceptible  to 
be  affected  by  ion  pairing. 

In  order  to  determine  precisely  the  B-  and  /9-coefficients 
in  the  JD  equation,  the  following  procedure  has  been 
employed.  A  plot  of  (rjv  —  1  )/C  against  C  forces  all  higher 
interaction  terms  to  be  absorbed  in  the  slope  [15,16]  and  the 
ordinate  at  C  =  0  gives  a  precise  value  of  B.  The  D-coeffi- 
cient  was  subsequently  determined  by  a  least  square  method. 
B-  and  D- values  obtained  by  this  method,  in  the  temperature 
range  25-55  °C,  are  given  in  Table  2. 

The  5-coefficients  vary  with  the  nature  of  the  anion  but 
not  with  the  temperature,  at  least  in  the  temperature  range 
investigated,  and  for  each  salt,  a  mean  value  has  been 
calculated.  As  indicated  by  the  S.D.  of  the  mean,  the 
observed  variations  are  in  the  order  of  the  experimental 
error.  At  the  opposite,  the  D-coefficients  decrease  regularly 
with  temperature  and  in  order  to  investigate  the  exact  nature 
of  the  temperature  dependency  of  D ,  we  have  plotted  D 
values  against  1  IT.  All  the  D  =f(l/T)  curves  reported  in 
Fig.  3  are  linear.  The  slopes  P  of  D  versus  1  IT  for  each  salt 
are  also  gathered  in  Table  2.  The  fact  that  D  is  more  sensitive 
to  thermal  agitation  than  B  is  not  surprising  since  B  is  mainly 
related  to  the  hydrodynamic  volume  (including  solvation 
effects,  refer  to  the  following  paragraph)  and  D  mostly  to 
ion-dipoles  interactions  and  long  range  coulombic  ion-ion 
interactions. 

The  BC  term  in  the  JD  equation  may  be  also  deduced 
from  the  hydrodynamic  equation  proposed  by  Einstein 
for  the  relative  viscosity  of  spherical  unsolvated  particles 

[17]: 

r]r  =  1  ±2.5  <P  +  o((P2) 


(7) 
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Fig.  2.  Variations  of  the  viscosity  of  LiBF4,  LiPF6,  LiAsF6  and  LiTFSI  solutions  in  BL  with  salt  concentration  at  25  °C. 


In  Eq.  (7),  (p  is  the  volume  fraction  of  the  particles  in  the 
fluid.  As  <P  can  be  expressed  as  a  function  of  the  concentra¬ 
tion  C  and  the  hydrodynamic  molar  volume  V)  of  the 
particles,  Eq.  (7)  may  be  rewritten  as: 

rjr  =  1  +  0.0025  U,C  +  o(C 2)  (8) 

when  Vi  is  in  cm  mol  and  terms  in  C  and  higher  are 
considered  as  negligible.  If  it  is  supposed  that  this  equation 
holds  at  the  molecular  scale  for  large  solute  molecules, 
identification  of  Eqs.  (1)  and  (8)  leads  to  Z?ca|c  =  0.0025  V,. 
Using  this  relation,  ones  can  deduce  5caic  from  the  solvated 
volume  of  the  electrolyte  V*.  Bcalc  values,  reported  in  Table  2, 
are  obtained  by  taking  r(Li+)soiv  =  0.346  nm  for  the  solva¬ 
ted  radius  of  the  Li  ion  [7,18]  and  r(anion)  =  0.227  nm 
(BF4  ),  0.254  nm  (PF6-),  0.259  nm  (AsF6-),  and  0.326  nm 
(TFSI-)  in  literature  [18].  The  agreement  between  experi¬ 
mental  and  calculated  values  of  B  is  good  which  means  that 
the  volume  effect  is  predominant,  at  least  for  large  anions 
and  strongly  solvated  cations. 


3.2.  Activation  energy  for  the  viscous  flow 

From  Eyring’s  theory  of  activated  processes  [19],  an 
expression  for  the  valuations  of  solutions  viscosity  with 
temperature  is: 


hNa  , 

(A 

\ 

\H+~ 

l 1  =  v  exp  ( 

^  m 

V  R  y 

Jexp 

RT 

In  Eq.  (9),  h  is  the  Planck’s  constant,  Vm  the  molar  volume  of 
solvent  (76.6  cm3  mol- 1  for  BL),  A the  activation  entropy 
and  A7U  the  activation  enthalpy,  generally  identified  to  the 
activation  energy  of  the  viscous  flow  Ean.  Non-associated 
solvents  and  non  glass-forming  ionic  and  molecular  liquids, 
usually  verify  this  equation. 

We  have  reported  in  Table  3  the  values  of  the  activation 
energy  and  entropy  for  the  viscous  flow  of  LiBF4,  LiPF6, 
LiAsF6,  and  LiTFSI  solutions  in  BL.  In  a  preceding  work  [7], 
a  linear  relation  between  the  activation  energy  for  the  viscous 


Table  2 


B  (AT 

')  and  D  (M  2)  coefficients  of  the  JD  equation  and  slope  (P/M 

2K)  of  the  D 

=/(  l/7)  relation  for  LiBF4 

h  LiPF6,  LiAsF6 

and  LiTFSI  solutions  in  BL 

T  (  C) 

LiBF4 

LiPF6 

LiAsF6 

LiTFSI 

B 

D 

B 

D 

B 

D 

B 

D 

25 

0.30 

0.75 

0.40 

1.10 

0.40 

1.15 

0.  50 

1.20 

35 

0.30 

0.70 

0.40 

0.95 

0.42 

1.00 

0.47 

1.10 

45 

0.30 

0.65 

0.38 

0.85 

0.42 

0.85 

0.46 

1.05 

55 

0.25 

0.55 

0.42 

0.68 

0.41 

0.75 

0.50 

0.90 

Anean 

0.288 

- 

0.400 

- 

0.413 

- 

0.483 

- 

±0.013 

±0.008 

±0.005 

±0.013 

#Calcd 

0.33 

- 

0.36 

- 

0.37 

- 

0.48 

- 

P 

630 

1330 

1320 

930 

a  Bca ic  =  0.0025  V[  (at  25  °C,  refer  to  text). 
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Fig.  3.  Variation  of  the  /7-coefficient  in  the  JD  relation  (Eq.  (1)  in  the  text)  with  1/71 


flow  Ea jj  and  the  salt  concentration  C  has  been  proposed: 

=  <„  +  V'n.EjC  (10) 

where  Ea  and  E^  respectively,  are  the  energy  of  activation 
for  the  pure  solvent  and  the  contribution  of  the  salt  (per 
mole  of  the  solute)  to  the  activation  energy  for  the  transport 
process. 


Table  3 

Activation  energy  Ea_n  (kJ  mol  ')  and  entropy  A (J  K  1  moP1)  for  the 
viscous  flow  of  LiC104,  LiBF4,  LiPF6,  LiTFSI  and  LiAsF6  solutions  in  BL 


C/M 

0 

0.1 

0.2 

0.5 

1.0 

1.5 

LiClO/ 

f 

a ,rj 

10.7 

11.1 

11.5 

12.6 

14.8 

16.3 

A  S* 

-12.2 

-11.4 

-10.8 

-9.3 

-5.5 

-4.7 

LiBF4 

f 

^ a.ri 

10.7 

13.1 

13.3 

14.4 

15.5 

17.2 

AS* 

-12.2 

-4.4 

-2.8 

-2.3 

-1.3 

-0.4 

LiPF6 

f 

10.7 

11.1 

12.0 

14.0 

17.7 

19.5 

AS* 

-12.2 

- 

-8.9 

-4.4 

1.9 

5.3 

LiTFSI 

L-,a,r\ 

10.7 

11.3 

11.7 

12.5 

14.5 

17.5 

A  S* 

-12.2 

-9.4 

-4.5 

-2.5 

0.2 

5.6 

LiAsF6 

f 

a.ri 

10.7 

12.0 

12.4 

12.7 

15.6 

17.3 

AS* 

-12.2 

1 

0° 

oo 

O' 

-8.1° 

-9.0 

-3.6 

-2.0 

a[7], 

bC  =  0.125  M. 
c  C  =  0.25  M. 


In  Fig.  4,  we  have  reported  the  variations  of  Ean  with  C  for 
the  different  salts  under  study.  Only  for  LiBF4,  the  ordinate 
at  the  origin  is  different  from  Ea  =  10.7  kJ  mol-1.  From 
the  slopes  of  the  lines,  E^f  values,  in  kJ  mol-1,  have  been 
infered:  38  (LiBF4),  8l’  (LiPF6),  60  (LiAsF6)  and  56 
(LiTFSI).  Since  Esaf  (LiC104)  =  49.3  kJ  mol-1  [7],  it 
appears  that  the  nature  of  the  anion  has  a  strong  influence 
on  Ef^,  BF4“  and  PF6  leading  respectively,  to  the  lowest 
and  the  largest  values.  One  may  notice,  that  the  most  the 
salt  is  expected  to  be  associated  in  ion-pairs  (especially  in 
the  case  of  small  anions  bearing  localized  charge),  the  most 
the  slope  and  hence  £^alt  is  low.  Ions  pairs  will  probably 
contribute  less  to  the  increase  in  viscosity  of  electrolyte 
solutions  than  free  ions. 

The  activation  entropy  A for  the  viscous  flow  can  be 
deduced  from  experimental  viscosities  using  Eq.  (9).  As 
shown  in  Table  3,  all  values  increase  when  the  salt  is  added 
and  at  C  =  1.5  M  follow  the  order: 

A5^(LiPF6)  ~  AS^  (LiTFSI)  >  AS^(LiBF4) 

>  AS^(LiAsF6)  >  AS^(LiC104) 

In  all  instances,  an  increase  of  AS/l  is  observed  when  the 
salt  is  added  to  the  pure  solvent.  This  means  that 
the  presence  of  the  salt  induced  a  decrease  in  order  in  the 
transition  state  relatively  to  the  initial  state.  As  it  is  con¬ 
sidered  that  the  presence  of  the  salt  at  high  concentration 
induces  the  formation  of  a  ionic  quasi-lattice  (vide  infra),  a 
decrease  in  entropy  and  consequently  a  sharp  increase  in 
order  in  the  initial  state  is  expected.  This  could  contribute 
to  explain  the  observed  increase  in  activation  entropy 
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Fig.  4.  Activation  energy  for  the  viscous  flow  (E.ul)  vs.  salt  concentration  (C)  for  LiBF4,  LiPF6,  LiAsF6  and  LiTFSI  solutions  in  BL. 

the  variations  of  A  for  LiBF4,  LiPF6,  LiAsF6  and  LiTFSI 
have  been  plotted  against  the  cube  root  of  the  salt  concen¬ 
tration  in  BL.  As  expected  from  Eq.  (3),  linear  plots  are 
obtained  for  all  salts  under  consideration. 

According  to  the  quasi-lattice  theory,  first  introduced  by 
GOSH  [10]  and  developed  further  by  Bahe  [20],  the  Debye 
length  kd  in  the  ionic  cloud  model  may  be  replaced,  when 
the  ionic  strength  is  raised  above  0.1  M  (as  for  these 
electrolytes),  by  /cL  the  average  distance  between  ions  of 
opposite  charge,  supposed  to  be  distributed  on  the  vertices  of 
an  expanded  lattice  [21]: 

kl  =  M(200(WaC)1/3  (11) 

In  Eq.  (11),  the  salt  concentration  C  is  expressed  in  mol  1  ’. 
In  a  similar  manner,  the  usual  DHO  conductivity  equation 
may  be  adapted  by  introducing  the  parameter  kl  instead  of 
kd  as  proposed  by  Smith  [22].  Then,  the  variations  of  the 
molar  conductivity  with  concentration  are  described  by 
Eq.  (3).  If  all  parameters  in  Eqs.  (3)-(6)  are  known,  S' 
for  a  given  salt-solvent  system  can  be  calculated  (S1^,)  and 
compared  to  the  experimental  value  (>S*xp).  S[.alc,  S'.xp,  A°' 
and  A°  values  are  gathered  in  Table  4.  S'.alc  has  been 
calculated  using  the  Madelung  constant  for  the  fee  lattice: 
M  =  1.74.  For  comparison,  values  obtained  for  LiC104  in 
BL  [7]  and  propylene  carbonate  (PC)  have  been  added  [23]. 
It  is  noticeable  that  values  are  close  to  SgXp  for  all  salts 
under  study,  as  previously  found  for  LiC104  in  BL  [7]. 
Nevertheless,  as  there  is  no  particular  reason  for  the  ions 
to  adopt  a  regular  fee  lattice,  the  constant  M  may  be  best 
regarded  as  an  adjustable  parameter.  Using  experimental 
data,  M  values  have  been  adjusted  to  experimental  data  and 
reported  in  Table  4.  All  values  fall  in  the  range  1.5-2. 4,  with 


(Fig.  5).  Nevertheless,  it  cannot  be  excluded  that  the  entropy 
of  the  transition  state  itself  is  affected  by  variations  in  ionic 
strength. 

3.3.  Conductivity 

3.3.1.  Variation  of  the  molar  conductivity 
with  concentration 

The  molar  conductivities  (A),  are  deduced  from  the 
experimental  conductivities  ( K)  by  A  =  K/C.  In  Fig.  6, 


Transition  state 


Initial  state 


AS* 


t  ransition  state 


Initial  state 


AS 


Transition  state 
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Fig.  5.  Possible  variations  of  the  entropy  of  activation  for  the  viscous  flow 
(A.U)  in  the  presence  of  a  dissolved  salt. 
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Fig.  6.  Plot  of  the  molar  conductivity  A  vs.  the  cube  root  of  the  concentration  C1/3  (C  expressed  in  mol  1  ')  for  LiBF4,  LiPF6,  LiAsF6  and  LiTFSI  solutions  in 
BL. 


a  mean  of  2.0,  not  to  different  from  1 .74,  the  rock-salt  value. 
A°'  and  A°  values  do  not  differ  significantly,  excepted  for 
LiBF4.  As  it  has  been  reported  that  the  dissociation  coeffi¬ 
cient  of  LiBF4  (1  M)  at  25  °C  in  BL  is  0.34  [24],  BL 
solutions  of  this  salt  may  involve  contact  ion-pairs,  leading 
to  a  possibly  underestimated  value  for  the  ordinate  in  the 
A  =/(C1/3)  plot. 

3.4.  Maximum  of  conductivity 

When  the  concentration  in  salt  is  raised,  the  number  of 
charge  earners  increases  but,  at  the  same  time,  the  viscosity 
increases  as  indicated  by  the  JD  relation.  The  consequence  is 
that  the  ionic  mobilities  inevitably  decrease.  The  competi¬ 
tion  between  the  increase  in  number  of  charge  carriers  and 
the  decrease  of  their  mobilities  lead  to  a  maximum  in  the 

Table  4 

Slopes  (S')  and  ordinates  ( A0 ')  in  the  correlation  A  =/(C1/3)  (Eq.  (3)  in  the 
text)  for  LiBF4,  LiPF6,  LiAsF6  and  LiTFSI  in  BL  and  LiTFSI  in  PC  at 
25  °C 


BL 

PC 

LiTFSI 

LiBF4 

LiPF6 

LiAsF6 

LiTFSI 

S' 

22.92 

36.55 

36.00 

26.47 

19.32 

s' 

calc 

23.54 

37.55 

36.34 

27.77 

11.65 

A0' 

30.29 

48.32 

46.77 

35.74 

24.69 

A0 

42.64a 

40.73“ 

39.91“ 

33.77“ 

23.80“ 

A°'/A° 

0.71 

1.18 

1.17 

1.06 

1.04 

M 

1.48 

2.40 

2.12 

1.86 

2.28 

A°  is  the  limiting  molar  conductivity  given  by  the  DHO  relation.  Eqs.  (4)- 
(6)  have  been  used  to  obtain  ,S(alc.  M  represents  the  value  of  the  optimized 
lattice  parameter. 
a  [18]. 


conductivity-concentration  relationship.  This  maximum  is 
usually  obtained  around  1  M  in  organic  electrolytes.  For 
evident  reasons,  the  determination  of  this  maximum  in 
conductivity  is  important  for  batteries  applications. 

The  maximum  conductivity  value  (Cmax)  is  easily 
deduced  from  Eq.  (3)  by  derivation  of  K  relatively  to  the 
concentration: 


Cmax  is  then  related  to  the  molar  conductivity  at  infinite 
dilution  A 0/  and  to  the  slope  of  the  A  =/(C1/3)  plot  which 
can  be  calculated  as  previously  indicated  (Eqs.  (4)-(6)). 

An  alternative  way  to  determine  Cmax  is  provided  by  the 
study  of  the  variations  of  the  Walden  product  W  =  A i]  with 
salt  concentration.  In  the  present  electrolyte  solutions,  W  is 
approximately  constant  as  indicated  by  the  Ln  //  =/(Ln  A) 
plot  reported  in  Fig.  7.  Straight  lines  are  obtained  whose 
slope  ranges  from  —0.8  to  —1  when  the  concentration  is 
varied  from  0.2  to  1.5  M.  Introducing  the  JD  equation  and  W 
in  the  expression  of  the  conductivity  of  electrolyte  solutions 
leads  to: 

WC  WC 

K  =  AC  = - =  (13) 

r\  (rf(l+BC  +  DC2))  V  ' 

Eq.  (13)  has  a  maximum  for: 

Cmax  =  JfTJl  (14) 

In  Table  5  are  reported,  the  concentration  at  maximum  of 
conductivity  obtained  from  experimental  data  and  values 
calculated  by  Eqs.  (12)  and  (14).  As  expected,  a  good 
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Fig.  7.  Ln  A  =/( Ln  ;/)  for  LiBF4,  LiPF6,  LiAsF6  and  LiTFSI  solutions  at  25  °C  in  BL. 


agreement  is  obtained  between  calculated  and  experimental 
values,  and  it  can  be  noted  that  Cmax  value  can  be  predicted 
using  only  A0'  as  experimental  entry  (approximated  by  A°  if 
necessary)  as  all  others  parameter  involved  in  S'  may  be 
calculated  (for  M  the  fee  value  is  used)  or  alternatively 
deduced  from  viscosity  measurements. 

As  the  quasi-lattice  model  disregards  the  formation  of 
non-conducting  contact  ion-pairs,  this  theory  cannot  be 
readily  extended  to  all  organic  electrolytes,  especially  those 
having  low  permittivity.  This  problem  will  be  considered  in 
a  following  paragraph. 

3.5.  Activation  energy  for  conductivity 

According  to  the  Eyring’s  theory,  viscosity  and  conduc¬ 
tivity  are  both  activated  transport  processes.  From  the  quasi¬ 
lattice  theory,  a  relation  between  the  activation  energy  for 


Table  5 


Salt  concentration  at  the  maximum  of  conductivity  (Cmax  in  M)  at  25  °C  in 
BL.  Eqs.  (12)  and  (14)  are  used  for  calculated  values 


Salts 

r 

'-^max 

cmaxa 

Cmaxb 

k  (mS/cm) 

LiBF4 

1.02 

0.90 

1.15 

7.33 

LiPFfi 

0.98 

0.90 

0.95 

11.63 

LiAsF6 

1.01 

0.90 

0.93 

10.48 

LiTFSI 

0.97 

0.90 

0.91 

9.21 

a  Given  by  Eq.  (12). 
b  Given  by  Eq.  (14). 


the  conductivity  (EaA)  and  the  salt  concentration  may  be 
inferred  [7]: 


Ea,A=E°aA+kC4/3 
where  k  is  given  by: 

2Aa/3e2VmM104 

(4ne0£rZ1/3) 


(15) 


(16) 


In  this  expression,  Z  is  the  number  of  ions  in  a  unit  cell  of  the 
cation  (or  anion)  sub-lattice.  EaA  versus  C473  has  been 
plotted  on  the  graphs  reported  in  Fig.  8.  As  previously  found 
for  the  BF-FiC104  system  [7],  the  variations  of  EaA  for 
FiBF4,  FiAsF6,  FiPF6  and  FiTFSI  follow  Eq.  (15).  Experi¬ 
mental  slopes  (k)  and  ordinate  (E®  ,)  values  obtained  by 
linear  regression  for  the  different  salts  are  reported  in 
Table  6.  The  ordinate  of  the  correlation  E®A,  which  repre¬ 
sents  the  infinitesimal  dilution  activation  energy  for  the 
conductivity,  reaches  a  value  which  is  close  to  the  activation 
energy  for  the  viscosity  of  pure  BL  (10.7  kJ  moE1).  At 
infinitesimal  dilution,  the  activation  energy  for  the  conduc¬ 
tivity  represents  the  interactions  of  one  mole  of  “dis¬ 
charged”  ions  (their  distance  is  considered  as  to  large  for 
even  long  range  interaction  to  occur).  Taking  M  =  1.74  and 
Z  =  4,  kcAc  given  by  Eq.  (14)  is  0.5  kJ  (mol/1)  4/3 .  There  is 
no  agreement  between  experimental  and  calculated  k  values. 
The  smallest  k  value  is  obtained  for  LiBF4.  This  lack  of 
agreement  indicates  that  others  interactions  than  coulombic 
forces  between  ions  of  opposite  charge  have  to  be  taken  into 
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account  to  obtain  a  quantitative  agreement  between  theory 
and  experience.  As  a  matter  of  fact,  the  ion-solvent  dipole 
interactions  energies  must  be  quite  high  in  media  containing 
strongly  dipolar  solvents  like  cyclic  carbonates  or  lactones. 

3.6.  Ion  association 


In  1926,  Bjerrum  [25]  introduced  the  notion  of  ion-pairs 
formed  by  the  following  chemical  equilibrium: 

A+  +  B~  (A+,B~) 


This  equilibrium  is  characterized  by  the  thermodynamical 
association  constant  K.d: 


(!  ~  «) 

(a2Cy±2) 


(17) 


In  Eq.  (17),  y±  is  the  mean  activity  coefficient  and  a  the  ion 
pair  dissociation  coefficient.  In  diluted  solutions,  Ka  is 
determined  using  the  DHO  treatment  for  a  and  the  Debye 
and  Hiickel  law  for  y±.  At  high  ionic  strength  (/  >0.05  M), 
the  use  of  the  preceding  equations  is  not  possible  and  the 
determination  of  both  a  and  y±  is  problematic.  For  these 
reasons,  the  ion  and  ion-pair  concentrations,  but  not  the 
activities,  are  best  determined  spectrocopically  using  IR  or 


Table  6 

Slopes  of  the  Elw\  =/(C4/3)  (kJ  moF1)  relationship  and  activation  energy 
at  infinite  dilution  (Eq.  (15)  in  the  text) 


Salts 

LiBF4 

LiPF6 

LiAsFf, 

LiTFSI 

k^M~m 

Ka 

1870 

10.8 

2120 

10.7 

2300 

11.0 

2690 

10.9 

Raman  spectroscopy.  Associated  species  like  contact  ion- 
pairs  (Li+. .  ,X  )  are  easily  detected  by  Raman  spectroscopy 
[26],  but  the  existence  of  solvent  separated  ions  pairs 
(Li+. .  .S  . .  .X-)  remains  questionable.  The  formation  of 
associated  species  depends  strongly  of  the  nature  of  the 
solvent  used.  As  an  example,  LiAsF6  remains  fully  disso¬ 
ciated  in  PC  [26]  but  in  dimethyl  carbonate  (DMC)  the 
following  species  have  been  detected  [26]:  free  Fi+  and 
AsFf)  ions,  contact  ion-pairs  [Fi+. .  ,AsF6  ],  solvated  ion- 
pairs  [Fi(DMC)4]+. .  .AsF6~  and  polymeric  ion-pairs 
[Fi+. .  ,AsF6_]v.  Even  in  high  permittivity  solvents,  contact 
ion-pairs  may  be  formed  when  the  salt  concentration  is 
increased.  This  is  in  particular  the  case  of  FiC104  in  PC  [27] . 

Ion  association  in  BE  has  been  investigated  for  LiPF6,  and 
LiC104  as  no  association  is  expected  for  the  TFSI~  anion 
which  has  the  largest  delocalized  structure.  In  fact,  it  has 
been  observed  by  Raman  spectroscopy  for  LiTFSI  in  BL 
[28]  that  unperturbated  BL  molecules  at  676  cm-1  coexist 
with  molecules  in  strong  interaction  with  Li+  ion  at 
690  cm  1 .  The  solvation  number  of  Li+  ions  found  by 
Raman  is  4  and  from  this  study,  it  results  that  only  free 
ions  and  possibility  solvent  separated  ions  pairs  may  exist  in 
LiTFSI/BL  solutions.  In  Fig.  9  is  reported  the  Raman  spectra 
of  solid  LiPF6  (as  powder)  and  LiPF6  in  BL  at  0.5  and  1.5  M. 
The  peaks  have  a  Gaussian  shape  and  the  area  (A)  under 
the  curves  is  a  linear  function  of  the  salt  concentration: 
A  =  0.0285C,  with  a  regression  coefficient  r  =  0.91.  There 
is  no  evidence  of  a  new  emission  band  which  could  be 
attributed  to  ion-pairing.  As  a  conclusion,  if  ions  pairs  are 
formed  at  high  concentration  in  salts,  these  would  be  solvent 
separated  or  solvent-solvent  separated  ion-pairs. 

Fig.  10  shows  the  Raman  spectra  of  BL  in  the  presence  of 
1  and  1.5  M  LiC104  in  the  range  915-955  cm-1  (after 
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Fig.  9.  Raman  spectra  of  solid  LiPF,,  (as  powder),  LiPF6  (0.5  M)  and  LiPF6  (1.5  M)  in  BL. 


subtracting  and  normalising  the  spectra  as  indicated  in  the 
experimental  part).  An  unperturbed  perchlorate  anion  has  a 
tetrahedral  structure  and  belongs  to  Td  symmetry  group.  The 
nine  vibrational  degrees  of  freedom  are  divided  into  four 
modes  of  vibration  [29]:  u i (A i )  =  931  cm-1,  112(^2)  = 
458  cm-1,  113(^2)  =  1100  cm-1,  tqC^)  =  624  cm-1.  All 
the  four  modes  are  Raman-active  and  F2  is  IR-active.  In 
LiC104  acetone  solutions,  the  following  species  have  been 
observed  [30]:  free  ions  CIO4  (933  cm-1),  contact  ions  pairs 
[Li+  . .  .CIO4  ]  (955  cm-1)  and  solvent  separated  ions  pairs 
(948  cm-1).  The  Raman  spectrum  of  LiC104  in  BL,  reported 
in  Fig.  10,  shows  the  presence  of  only  two  emission  bands 


which  have  been  attributed  to  free  ions  and  ion-pairs.  The 
exact  nature  of  the  ion-pairs  is  not  known.  Using  the  emission 
at  the  preceding  wavelength,  the  dissociation  coefficient  a  of 
LiC104  in  BL  has  been  evaluated  as  a  function  of  the  salt 
concentration.  The  a  values,  reported  in  Table  7,  exhibit  little 
variations  when  the  salt  concentration  of  the  is  raised  from  0. 1 
to  1.5  M  as  the  precision  of  the  measurements  is  of  the  order  of 
±0.05.  As  a  consequence,  the  cube  root  law,  Eq.  (3),  remains 
applicable  and  only  the  slope  of  the  A  =f(C  )  relation  will 
be  affected  by  the  presence  of  ion-pairs. 

The  formation  of  loose  ion-pairs  (solvent  separated  or 
solvent-solvent  separated  ion-pairs)  is  highly  probable  even 


Fig.  10.  Raman  spectra  of  BL  without  salt  and  with  added  LiC104  (1  M  and  1.5  M)  in  the  915-955  cm  1  range. 
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Table  7 

Dissociation  degrees  of  LiC104  in  BL  at  25  °C  determined  by  Raman 
spectroscopy 


C/M 

a  (LiC104) 

0.10 

0.54 

0.20 

0.46 

0.50 

0.67 

0.75 

0.64 

1.00 

0.65 

1.50 

0.55 

in  high  permittivity  organic  electrolytes,  but  their  detection 
by  Raman  spectrometry  cannot  be  always  achieved,  as  the 
emission  bands  of  anions  are  only  sligthly  affected  by  ion 
pairing  in  this  case.  In  any  case,  their  presence  is  not  in 
conflict  with  the  quasi-lattice  structure  of  the  electrolyte  as 
these  species  are  only  temporarily  formed  and  can  be 
considered  as  defects  in  the  quasi-lattice  structure  of  the 
electrolyte  solution. 

4.  Conclusion 

The  aim  of  this  work  is  to  highlight  the  transport  proper¬ 
ties,  viscosity  and  conductivity,  of  concentrated  organic 
electrolytes.  The  duality  between  A  and  rj,  found  in  the 
Walden  product  W  for  example,  find  an  easy  interpretation 
in  terms  of  the  transition  state  theory  when  the  activation 
energies  for  conductivity  and  viscosity  have  similar  values 
[31].  The  quasi-lattice  model,  widely  used  in  fused  salt 
studies,  can  be  successfully  extended  to  concentrated  organic 
electrolytes  involving  strongly  dipolar  aprotic  solvents. 

The  coulombic  interactions  between  neighbouring  ions  on 
which  this  model  is  based,  are  able  to  explain  the  variation  of 
the  molar  conductivity  A  and  its  activation  energy  Ea  A  when 
the  concentration  of  the  salt  is  varied.  The  same  forces  are 
at  the  origin  of  the  increase  in  viscosity  //  of  the  electrolyte 
which  is  well  described  by  the  JD  semi-empirical  law.  The 
bell  shape,  usually  observed  for  the  K  =/(C)  plot  in  dipolar 
aprotic  solvents  is  well  in  accordance  with  the  cube  root  law 
for  A.  The  predictive  character  of  the  model  is  potentially 
interesting  as  the  only  adjustable  parameter  M,  the  lattice 
parameter  has  been  found  to  vary  only  moderately  from  the 
fee  value.  Combination  of  the  Walden  product  and  of  the  JD 
equation  shows  that  the  concentration  in  salt  at  the  max¬ 
imum  of  conductivity  Cmax  is  strongly  related  to  the  D- 
coefficient  and,  hence  to  ion-ion  interactions. 

The  association  of  ions  in  pairs,  as  observed  by  Raman 
spectroscopy  for  LiC104  in  BL,  does  not  invalid  the  use  of 
the  cube  root  law  as  it  has  been  found  that  the  dissociation 
coefficient  remains  constant  in  the  concentration  range  0.5- 
1  M.  Others  associated  species  than  contact  ion-pairs  may  be 
considered  as  defects  in  the  quasi-lattice  structure  of  the 
electrolyte  solution. 
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